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ABSTRACT 

To search for circumstcllar disks around evolved stars, we targeted roughly 100 
DA white dwarfs from the Palomar Green survey with the Peters Automated Infrared 
Imaging Telescope (PAIRITEL). Here we report the discovery of a debris disk around 
one of these targets, the pulsating white dwarf PG 1541+651 (KX Draconis, hereafter 
PG1541). We detect a significant flux excess around PG1541 in the iGband. Follow- 
up near-infrared spectroscopic observations obtained at the NASA Infrared Telescope 
Facility (IRTF) and photometric observations with the warm Spitzer Space Telescope 
confirm the presence of a warm debris disk within 0.13-0.36 Rq (11-32x the stellar 
radius) at an inclination angle of 60°. At T e g = 11880 K, PG1541 is almost a twin 
of the DAV white dwarf G29-38, which also hosts a debris disk. All previously known 
dusty white dwarfs are of the DAZ/DBZ spectral type due to accretion of metals from 
the disk. High-resolution optical spectroscopy is needed to search for metal absorption 
lines in PG1541 and to constrain the accretion rate from the disk. PG1541 is only 55 
pc away from the Sun and the discovery of its disk in our survey demonstrates that 
our knowledge of the nearby dusty white dwarf population is far from complete. 

Key words: white dwarfs — stars: individual (PG 1541+651, KX Dra) — infrared: 
stars — infrared: planetary systems 



1 INTRODUCTION 



Eighteen years a fter the discovery of the first dusty white 



Bulk Earth JZuckerman et al.ll2007l: iKlein et al.ll20ld . l201ll : 



iDufour et alJl201Ch 



dwarf (G2 9-28. IZuckerman fc Becklinl 1 19871 ). iKilic et all 
l|2005l ) and lBecklin et all l|2005l ) identified the second white 
dwarf with a circumstellar debris disk, GD 362. The lat- 
ter w as the most metal-ric h white dwarf known at the 
time (Giannin as et all 120041 ') . The discovery of circumstel- 
lar disks around these two DAZ stars led to a paradigm 
shift in our understanding of the source of the metals in 
white dwarf photospheres (|Kilic fc Red ficld 2007j) and also 
fueled further searches for debris disks around metal-rich 
white dwarf s. The discovery of the third (|Kilic et al.ll2006h 
and fourth (|von Hippel et all 2007) dusty white dwarfs fol- 
lowed shortly. Thanks to the Spitzer Space Telescope, we 
now know of >20 white dwarfs with circumstellar disks 
(seelFarihi et al.ll20ld,l201ll:lDufour et al.|[201ol:lDebes ; et"afl 
20111: iMelis et al.ll201ll : IXu fc Jurall201ll : iGirven et ai]|201ll : 



Steele et alj|20l3 . and references therein). A detailed photo- 



spheric abundance analysis of these white dwarfs show that 
the accreted metals are likely to originate from tidal disrup- 
tion of minor planets with compositions similar to that of 



All of the previously known dusty white dwarfs show 
metal-absorption features in their optical spectra (DAZ or 
DBZ spectral type), and they were specifically targeted (ex- 
cept G29-38) for infrared searches because they are metal- 
rich. The best way to identify new dusty white dwarfs 
would be to observe more DAZ/DBZ white dwarfs in the 
infrared. However, the classification of a DA as a metal- 
rich star requires the detection of weak metal lines with 
equivalent widths of ~ 50 mA. This can only be achieved 
with high resolution spectroscopy at 10m class telescopes 
|Zuckerman et a l. 2003; Koester et al.l l2005h A cheaper way 
to find more dusty white dwarfs is to perform a near-infrared 
search around a large sample of white dwarfs without prior 
knowledge of metal absorption features in the targeted stars. 
With high-quality optical and near-infrared photometry, a 
warm debris disk sign ature can be recognized in the K- 
band |Kilic et alj|2006h . For example, a K-band only search 
would detect the disks around 70% of the known dusty white 
dwarfs. 
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To identify new dusty white dwarfs and to constrain 
the frequency of dusty white dwarfs among the DA popu- 
lation, we performed a near- infrared JHK survey of «100 
DA white dwarfs from the Palomar Gree n (PG) survey usin g 
the robotic 1 . 3m te lescope PAIRITEL dBloom et alj|2006h . 
iLiebert et al.l (|2005l ) performed a detailed model atmosphere 
analysis of the DA white dwarfs in the PG survey using op- 
tical spectroscopy and provide temperature, surface gravity, 
mass, and age estimates. Their spectroscopy does not have 
enough resolution to detect the metal-rich DAZs. We select 
fslOO apparently single white dwarfs with T cff = 9000-22,000 
K, where we are most efficient in finding the disks using 
near-infrared data. 

Here we report the discovery of a circumstellar debris 
disk around one of our targets, the DA Variable (DAV) white 
dwarf PG1541. The rest of our sample will be presented in 
an upcoming publication. In Section 2 we describe our in- 
frared photometric and spectroscopic observations. In Sec- 
tion 3 we present the spectral energy distribution of PG1541 
and constrain the physical parameters of the circumstellar 
disk. We discuss the population of dusty white dwarfs in the 
PG survey in Section 4. 




Figure 1. The SDSS g-band and PAIRITEL J (top panel), H 
and K B (middle panel), and IRAC 3.6 and 4.5/^m (bottom panel) 
images of PG1541. 



2 OBSERVATIONS 

2.1 Near-Infrared Photometry 

PG1541 is a 510 Myr old ZZ Ceti variable white dwarf with 
T cff = 11880 K log. q = 8.2, M = 73Mb, R = 0.011-R S 
and d — 55 pc (Gi anninas et al]|201ll ). It was originally an 
A type star with an i nitial mass of 3 M l|Kalirai et al.ll200Sl ; 
IWilliams et~ai]|2009h . It has an apparent magnitude of g = 
15.6 in the Sloan Digital Sky Survey (SPSS ). It is detected 
in th e Two Micron All Sky Survey (2MASS. ISkrutskie et all 
l2006h with J = 15.604 ± 0.062, if = 15.912 ± 0.171, and 
K s = 15.429 ± 0.175 mag. Due to inaccurate H- and K B - 
band photometry, it is impossible to search for an infrared 
excess around this object using the 2MASS data alone. 
We obtained simultaneous JHK S imaging of PG1541 us- 
ing the PAIRITEL on UT 2010 February 15. PAIRITEL 
is the old 2MASS telescope and it uses the same cam- 
era and the filter set as 2MASS. We use the PAIRITEL 
Pipeline version 3.3 processed images and standard IRAF 
DAOPHOT routines to perform aperture photometry on 
every 2MASS source detected in the images. We use the 
2MASS stars to calibrate the photometry for PG1541. We 
measure J = 15.625 ± 0.013, H = 15.590 ± 0.025, and 
K B = 15.450 ± 0.048 mag in the 2MASS magnitude sys- 
tem. Our J- and Jf s -band photometry is consistent with 
the 2MASS measurements within the errors. Thanks to the 
smaller errors in the PAIRITEL photometry, a significant 
flux excess is revealed in the K s -ba,nd. 



Table 1. Optical and Infrared Photometry of PG1541 



Filter 


A eff (/Mn) 




FUV 


0.15 


574.4 ± 25.3 


NUV 


0.23 


1330.3 ± 25.7 


u 


0.36 


1567.4 ± 29.4 


9 


0.47 


2186.9 ± 41.5 


r 


0.62 


1971.6 ± 32.4 


i 


0.75 


1677.6 ± 22.4 


z 


0.89 


1273.1 ± 16.9 


J 


1.24 


896.4 ± 19.0 


H 


1.66 


594.7 ± 18.0 


K s 


2.16 


440.5 ± 21.2 


IRAC1 


3.56 


435.4 ± 13.5 


IRAC2 


4.51 


483.9 ± 15.0 



resolving power of 90-210 over the 0.7-2.5 /jm range. The 
observations are taken in two different positions on the slit 
separated by 10". The total exposure time for PG1541 is 32 
min. We use internal calibration lamps (a 0.1W incandescent 
lamp and an Argon lamp) for flat-fielding and wavelength 
calibration, respectively. To correct for telluric features and 
flux calibrate the spectra, we use the observations of the 
nearby bright A0V star HD 143187 at an airmass similar 
to the PG1541 observati ons. We use the ID L-based pack- 
age Spextool version 3.4 l|Cushing et al.|[2004l ) to reduce the 
data. 



2.2 Near-Infrared Spectroscopy 

We obtained low resolution near infrared spectra of PG 1541 
on 2011 April 13 using the 3m NASA Infrared Telescope 
Facility (IRTF) equipped with the 0.8-5.4 Mi cron Medium- 
Reso lution Spectrograph and Imager (SpeX, iRavner et al.l 
l2003h . The observing setu p and procedures are similar to 
those of iKilic etafl j2008l ). We use a 0.5" slit to obtain a 



2.3 Spitzer Photometry 

We used the war m Spitzer equipped with the InfraRed Array 
Camera (IRAC. iFazio et al.l 12004) to obtain infrared pho- 
tometry of PG1541 on UT 2011 January 6 as part of the 
program 70023. We obtained 3.6 and 4.5/im images with 
integration times of 30 seconds for nine dither positions. 
Figure 1 shows the optical and infrared images of the field 
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around PG 1541. These images reveal no contaminating 
sources around PG1541. We use the IDL astrolib packages 
to perform aperture photometry on the individual corrected 
Basic Calibrated Data frames from the S18.18.0 pipeline re- 
duction. We get similar results using 2, 3, and 5 pixel aper- 
tures, but we quote the results using the smallest aperture 
since it has the smallest errors. 

Following the IRAC calibration procedures, we correct 
for the location of the source in the array before averag- 
ing the fluxes of each of the dithered frames at each wave- 
length. We also correct the Channel 1 (3.6/im) photome- 
try for the pixel-phase-dependence. We estimate the pho- 
tometric error bars from the observed scatter in the nine 
images corresponding to the dither positions. We add the 
3% ab s olute calibra t ion e rror in quadrature dReach et al.l 
l2005al ). lReach" et al. | <|2009h demonstrate that the color cor- 
rections for dusty white dwarfs like G29-38 are small (0.4- 
0.5%) for channels 1 and 2. We ignore these corrections 
for PG1541. We measure F„ = 435.4 ± 13.5/^Jy in Chan- 
nel 1 and 483.9 ± 15.0/iJy in Channel 2. Table 1 presents 
the G alaxy Evolution Explorer (GALEX, iMorrissev et al.l 
120071 ) ultraviolet, SDSS optical, and PAIRITEL / Spitzer in- 
fr ared photometry of PG 1541. We use the corrections given 
in lEisenstein et all (|2006t ) to convert the SDSS photometry 
to the AB magnitude system. 



3 RESULTS 

Figure 2 presents the flux calibrated spectrum of PG1541 
along with the optical photometry from the SDSS and the 
PAIRITEL near-infrared photometry. The observed spec- 
trum is shown as a jagged line, and the expected near- 
infra red photospheric flux from the star ([Bergeron et al.l 
Il995h is shown as a solid line. The model spectrum for 
PG1541 is ca lculated by one of us (P .P.) using the param- 
eters given in lGianninas et all (|201ll ). To match the resolu- 
tion of the IRTF spectrum, the model white dwarf spectrum 
is shown at 100 A resolution. The dotted line shows the tel- 
luric spectrum obtained from the reference A0V star obser- 
vations. The telluric correction for PG1541 is not perfect in 
the regions around 1.4 and 1.9 /im. 

The IRTF spectrum matches the observed photometry 
in the 0.7-2.5/im region relatively well and it confirms the 
slight //-band excess and more significant if s -band excess in 
the PAIRITEL photometry. The shape of the infrared excess 
is similar to the excess seen in the other dusty white dwarfs; 
there are no obvious emission or absorption features, and it 
is rising toward longer wavelengths. 

Figure 3 shows the spectral energy distribution of 
PG1541 in the 0.1-5,um range, including the GALEX, SDSS, 
PAIRITEL, and Spitzer photometry. PG1541 is relatively 
nearb y and the line of sig ht reddening is only E(B — V) = 
0.03 (Schl egel et al.l [l998h ; its spectral energy distribution 
does not show any evidence of reddening. 

PG1541 is brighter at 4.5^im than it is at 2[im. Since 
the stellar contribution falls like a power law in this wave- 
length range, the contribution from the cool component is 
still rising at 5/im. This contribution is almost identical to 
the excess flux seen aro und G29-38 in the same wavelength 
region (see Figure 3 in iReach et al.l I2009T ) . A brown dwarf 
companion cannot explain the observed photometry in the 
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Figure 2. Flux calibrated spectrum of PG1541 (jagged line) and 
the telluric and instrumental spectrum derived from the reference 
A0V star (dotted line). The expected flux from the stellar pho- 
tosphere is shown as a solid line. SDSS optical photometry and 
PAIRITEL near-infrared photometry are shown as filled points. 



2-5 (jm range l|Leggett et alj|20ld ). In addition, spectral fea- 
tures from such a companion would be visible in our IRTF 
spectrum. Hence, the infrared excess around PG1541 can 
only be explained by a circumstellar disk. 

We fit the infrared excess ar ound PG151 4 using the op- 
tically thick flat-disk models of I Jural ((2003) . Using the ef- 
fective temperature, radius, and distance estimates for the 
white dwarf, we created a grid of disk models with inner 
temperatures 1000-1800 K, outer temperatures 300-1300 K 
(in steps of 100 K and with the condition that Tin > T ou t) 
and inclination angles 0° — 90° (in steps of 10°). The model 
that best fits the infrared photometry has Tin = 1300 K, 
Tout = 600 K, and an inclination angle of i — 60°. To es- 
timate the uncertainties for these parameters, we perfom a 
Monte Carlo analysis where we replace the observed pho- 
tometry / with / + gSf, where Sf is the error in flux and g 
is a Gaussian deviate with zero mean and unit variance. For 
each of the 10,000 sets of modified photometry, we repeat 
our analysis and derive new best-fit parameters for the disk. 
We adopt the interquartile range as the uncertainty. This 
analysis yields T in = 1300 ± 100 K, T out = 600 ± 100 K, 
and i = 60 ± 10 degrees. The disk is within 0.13-0.36 Rq 
of the white dwarf; this i s within the Roche rad ius for tidal 
disruption of an asteroid (|von Hippel et al1l2007h . These pa- 
rameters are similar to the parameters for the other dusty 
white dwarf s studied w ith the same flat-disk model, includ- 
ing G29-38 (|Jurall2003l ). 

The total mass of the disk cannot be constrained for 
an opaque ring. Based on the choice of optically thin or 
thick mo dels, the di s k around G29-38 ho lds 10 19 - 10 24 g of 
material (|Jurall2003l : IReach et alj|2005bh . A similar amount 
of material may exist in the disk around PG1541. 
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Figure 3. Spectral energy distribution of PG1541 from ultravio- 
let to infrared wavelengths. The GALEX, SDSS, PAIRITEL, and 
IRAC photometry as well as our IRTF spectrum are shown. The 
solid line represents the contribution from the stellar photosphere 
and the open circles represent the average over the filter band- 
passes. The dotted line shows the contribution from the debris 
disk. 



4 DISCUSSION 

Along with GD 133, WD 1150-153, and G29-38, PG1541 
is the fourth DAV known to host a circumstellar debris 
disk. Every DA is expected to go through the ZZ Ceti in- 
stability strip during its evolution (see the discussion in 
iGianninas et al.ll201ll ). Hence, we expect that some of the 
DAZs with disks will also be DAVs; there is nothing special 
about finding 4 dusty DAV white dwarfs. However, optical 
and infrared observations of the pulsations in the first known 
dusty white dwarf, G29-38, were important for identifying 
the source of the infrared e xcess as a debris disk rather than 
a brown dwarf companion l|Graham et al.lll99fl) . 

The main pulsation mode in PG1541 ha s an amplitude 
of 4.5% in the optical wit h a period of 689 s (|Vauclair et alj 
l2000l ). lReach et al.l (|2009f i detect 4% fluctuations in G29-38 
at 3.6/im, whereas the main pulsatio n mode in the optical 
has a n amplitude larger than 20% l|McGraw fc Robinson! 
1 19751 ). The infrared variations are likely due to the dust 
grains going through temperature fluctuations with the stel- 
lar pulsations. Assuming that the disk around PG1541 is 
similar to G29-38's disk, we do not expect to see any sig- 
nificant (>1%) variations in our IRAC data for PG1541. 
Hence our disk model is not likely to be effected by stellar 
variability. 

Previous Spitzer su rveys find a freque n cy of 1-3% 
for dusty white dwarf s dMullallv et al 1 120071 ; iFarihi et alj 
20091; iKilic et al.ll2009r i. Out of the 348 DAs analyzed by 



Liebert et all (|2005T ). there are now six PG white dwarfs 
known to host debris disks: PG 0843+517, 1015+161, 
1116+026, 1457-086, 1541+651, and 2326+049 (G29-38). 



Therefore, the frequency of debris disks in the PG sample 
is at least 1.7%. However, not all of the DA white dwarfs 
in the PG survey have been observed in the infrared, and 
a more accurate analysis will have to wait until our entire 
sample is presented in an upcoming publication. 

Th e ongoing Wide-fi eld Infrared Survey Explorer 
(WISE, IWright et al.ll2010Ti mission should greatly increase 
the sample of known dusty white dwarfs in the solar neigh- 
borhood. The WISE InfraRed Excesses around Degenerates 
(WIRED) survey has already identified a new d usty white 
dwarf at 55 pc, GALEX1931 l|Vennes et al.ll2010ri . and it is 
predicted to detect disks around G29- 38-like objects within 
~100 pc of the Sun (|Debes et al.ll20l"lh . However, due to its 
relatively large beam size (> 6.1"), background contamina- 
tion may be a problem for WISE observations of targets in 
crowded fields (jMelis et alj|201ll ). PG1541 is bright enough 
to be detected in the WISE 3.4 and 4.6 /jm bands, and possi- 
bly in longer wavelength data if it displays a silicate emission 
feature. 



5 CONCLUSIONS 

We have detected an infrared excess at 2-5 /im around 
the 510 Myr old (cooling age) pulsating white dwarf PG 
1541+651 and determined it to be a narrow dust ring within 
the tidal disruption radius of the white dwarf. Near-infrared 
spectroscopic data show a significant contribution from a 
featureless spectrum rising to the red, confirming our dust 
disk interpretation. PG1541 is most likely a DAZ, however 
high-resolution optical spectroscopy is required to identify 
weak metal absorption features. 

We detected the disk around PG1541 in our survey of 
~100 DA white dwarfs from the PG survey. Our remaining 
targets will be presented in a future publication, which will 
also address the frequency of disks in the PG survey. 
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